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! Abstract 

The tight couphng expansion, appropriately generahzed to include large-scale magnetic 
^ ', fields, allows the estimate of the brightness perturbations of CMB anisotropics for typical 
wavelengths that are larger than the Hubble radius after matter-radiation equality. After 
discussing the basic features of the the pre-decoupling initial conditions in the presence of 
fully inhomogeneous magnetic fields, the tight coupling expansion is studied both analyti- 
cally and numerically. From the requirement that the amplitudes and phases of Sakharov 
oscillations are (predominantly) adiabatic and from the inferred value of the plateau in the 
temperature autocorrelation, the effects of the magnetized contribution can be systematically 
investigated and constrained. 
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1 Introduction 



A distinctive signature of pre-decoupling (adiabatic) initial conditions is that the temperature 
autocorrelation has the first (Doppler) peak for — 220 [1, 2, 3]. For the same set of initial 
conditions, the cross-correlation power spectrum between temperature and polarization will 
have the first (anticorrelation) peak for 4 — 3 -^(1/4 < id [1]- 

Defining as k the comoving wave-number and as Csb the sound speed of the baryon- 
photon system, the large-scale contribution to the temperature autocorrelation oscillates as 
a cos (A;CsbTdcc)j while the cross-correlation oscillates as sin (2A;CsbTdcc)- The first (compres- 
sional) peak of the temperature autocorrelation will then correspond to kCshTdec ~ (i.e. 
£d — 220), while the first peak of the cross-correlation will arise for kcshTdec ~ 37r/4 (i.e., 
approximately, ic ^ 3^d/4). This result can be obtained analytically by working to first- 
order in the tight-coupling expansion [4, 5, 6]. In this framework, the key assumption is 
that the initial conditions of the Einstein-Boltzmann hierarchy are characterized, prior to 
matter-radiation equality, by a solution that is, at least predominantly, adiabatic. A corol- 
lary to the previous statements is that the peak in the cross-correlation occurs for typical 
wavelengths that are larger than the Hubble radius at the the time of photon decoupling. 

In a previous paper [7] it has been shown that large-scale magnetic fields, present prior 
to equahty affect the CMB initial conditions. If only the adiabatic mode is present, the fully 
inhomogeneous magnetic fields contribute to the Sachs- Wolfe plateau as a subdominant 
non-adiabatic component that may be either correlated or uncorrelated with the adiabatic 
mode. If the pre-decoupling initial conditions are given by a mixture of adiabatic and non- 
adiabatic modes, the magnetized contribution may also mix with the isocurvature modes so 
that the total number of parameters defining the initial conditions increases. It is desirable 
to complement and extend the analysis of [7] by scrutinizing the impact of pre-decoupling 
magnetic fields on intermediate scales. Within the tight coupling expansion, appropriately 
generalized to include the effects of fully inhomogeneous magnetic fields, the fiuctuations in 
the brightness induced by the (adiabatic and non-adiabatic) scalar modes can be estimated. 

The effects of magnetic fields on the scalar brightness perturbations are typically ne- 
glected. To give an example, let us consider how the impact of large-scale magnetic fields 
is usually computed, for instance, in the Faraday effect. This case is interesting, since, the 
rotation of the CMB polarization plane indeed presupposes that the magnetic fields do not 
affect the process of formation of polarization. The influence of large-scale magnetic flelds on 
CMB polarization can be separated, for illustration, into two physically distinct effects. Prior 
to decoupling, gravitating magnetic fields modify the initial conditions of CMB anisotropics 
and affect the evolution equations of the photon-baryon system. Therefore, the induced de- 
gree of polarization will be sensitive to the presence of magnetic fields whose role will be, at 
this stage, to modify the features of the adiabatic initial conditions. This is the first effect. 
Once the polarization is formed, the polarization plane of the CMB may be Faraday rotated 
[8, 9]. From this second effect interesting bounds on the magnetic field intensity are usually 
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derived. 

This philosophy motivated various independent studies on the possible Faraday rotation 
of the CMB polarization degree [10, 11, 12, 13]. In some cases the magnetic field has been 
assumed just uniform [10], but there are also calculations analyzing the situation where the 
magnetic field is fully inhomogeneous [11, 12]. In spite of the relevant technical differences, 
the common assumption of the mentioned investigations is that the magnetic field is negligi- 
ble in the process of formation of CMB polarization. The latter assumption can be rephrased 
by saying that the initial conditions of the Einstein-Boltzmann hierarchy are set in such a 
way that the magnetic fields are absent. In the absence of large-scale magnetic fields the only 
source of polarization arises to first-order in the tight coupling expansion. The magnitude of 
the obtained contribution is just a consequence of the adiabaticity of the solution. If large- 
scale magnetic fields are consistently introduced before matter-radiation equality, the usual 
adiabatic (as well as non-adiabatic) initial conditions of the Einstein-Boltzmann hierarchy 
are modified and, consequently, also the degree of polarization of the CMB is affected. 

This example motivates a systematic analysis of the impact of magnetized initial con- 
ditions on the brightness perturbations of the radiation field. As anticipated above, this 
calculation will be conducted within the tight-coupling approximation which is known to be 
rather effective in obtaining the brightness perturbations for large and intermediate scales. 
This analysis has never been done before, to the best of our knowledge. A related motivation 
is that, in the absence of magnetic fields, the tight coupling approximation is employed in 
numerical Boltzmann codes at early times to avoid the integration of stiff Euler equations. 
The results of the present investigation allow to use the same strategy also in the presence 
of fully inhomogenous magnetic fields. 

Several aspects of CMB physics are germane to large-scale magnetization (see [14, 15] for 
two topical review articles). The astrophysical evidence for large-scale magnetization can be 
summarized as follows. Large-scale magnetic fields are an important element of the physics 
of the interstellar medium and have been measured, through various techniques in galaxies 
[16, 17] and in rich clusters [18, 19, 20, 21]. The determination of large-scale magnetic fields 
associated with loosely gravitationally bound systems (like superclusters) is still debatable 
[22, 23, 24], but, nonetheless, extremely intriguing. As far as galactic magnetic fields are 
concerned, there is the hope that in the future all sky survey of Faraday rotation (achievable 
through the Square Kilometer Array [SKA] [17, 18]) will allow a Faraday "tomography" 
of the Milky Way addressing the dark corners (number of reversals, pitch angle) of the 
morphological structure of the closest large-scale magnetic field we observe through different 
surveys of pulsars such as the Parkes and the southern galactic plane surveys (see [25, 26, 27] 
and references therein). There is the hope, in the present framework, that the present-day 
coherent field in spiral galaxies may keep the "memory" of the structure of the initial seed 
field (see [16] for a review of this long-standing problem). The morphological features of the 
magnetic fields profiles in clusters are less understood than in the case of galaxies. However, 
attempts have been made in constraining the radial profiles of magnetic fields in clusters 



2 



[19, 20]. Determinations of intracluster magnetic fields through statistical samples of rich 
clusters [21] seem to suggest that (ordered?) intra-cluster magnetic fields may be as large as 
the /xG (see, for instance, [28] for a review on the features of our magnetized Universe). 

Compressional amplification (taking place during the gravitational collapse of the proto- 
galaxy) allows to connect the observed field to a protogalactic field, present prior to gravita- 
tional collapse, between 0.1 nG and 0.01 nG. It is plausible, within the dynamo hypothesis, 
that the protogalactic fields could even be much smaller than the nG. In recent years, a lot 
of progress has been made both in the context of small- and large-scale dynamos i.e. the 
process by which the kinetic energy of the protogalactic plasma is transferred, by differential 
rotation, into magnetic energy. This progress [29] (sec also [30] for an introduction to astro- 
physical dynamos) has been driven both by the higher resolution of numerical simulations 
and by the improvement in the understanding of the largest magnetized system that is close 
to us, i.e. the sun. In contrast with the previous lore, it is now clear that the dynamo action 
demands a change in the topology of the fiux lines. As a consequence, large-scale dynamos 
should also produce small-scale helical fields that quench (i.e. prematurely saturate) the a 
effect [29] (see also [31]). The morphology and strength of magnetic fields in clusters may 
be related to the way the dynamo effect saturates. 

The simplest way to understand why large-scale magnetic fields may be relevant for CMB 
physics is to think of a protogalactic magnetic field in the nG range. If this field is naively 
blue-shifted at the decoupling epoch, i.e. at a redshift ^dec — HOO its strength could be as 
large as the mG, i.e. roughly six orders of magnitude larger. In the past ten years, indeed 
a lot of work has been done with the purpose of constraining large-scale magnetic fields 
using CMB physics. Historically the first type of configurations used in this type of exercise 
have been weakly uniform magnetic fields that would partially break the isotropy of the 
background geometry [15]. It was than realized that this assumption induces a correlation 
of the a^-i,m and a^+i,™ multipole coefficients arising in the angular power spectrum (see 
[32, 33, 34] and references therein). From this observation, uniform magnetic fields can be 
constrained. 

Another (complementary) possibility is that magnetic fields are fully inhomogeneous. In 
this case the magnetic fields do not break the spatial isotropy of the background but they 
can affect virtually all observables relevant in the theory of CMB anisotropics. In particular, 
in a series of papers. Barrow and Subramanian [35, 36, 37] pointed out possible effects of 
magnetic fields on the induced vector CMB anisotropics as well as on the polarization power 
spectra induced by the same modes. Within a slightly different approach the vector modes 
(as well as the tensor modes) have been discussed in [38]. More recently [39] a full numerical 
analysis of the vector and of the tensor modes induced by fully inhomogeneous magnetic 
fields has been presented. In [40] the signatures of magnetic helicity have been studied 
always in the perspective of the vector modes. Finally, in a series of papers, the evolution 
of weakly inhomogeneous magnetic fields has been discussed within the covariant formalism 
and in the presence of a uniform magnetic background [41]. 
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In spite of the fact that vector modes are probably one of the characteristic signatures 
of large-scale magnetic fields, scalar modes are also induced [7, 42]. This statement can be 
understood by noticing that the magnetic energy density, the magnetic anisotropic stress and 
the magnetic pressure contribute, respectively, to the Hamiltonian and momentum constraint 
and to the dynamical evolution of the curvature perturbations. Furthermore, within the 
magnetohydrodynamical approximation, the divergence of the Lorentz force contributes to 
the evolution equation of the divergence-full part of the baryon peculiar velocity. 

What is missing, at the moment, is the analysis of the brightness perturbations induced 
by magnetized initial conditions. As anticipated, one of the purposes of the present paper 
will exactly be to bridge fill this gap and to allow an approximate evaluation of the brightness 
perturbations, also at intermediate scales, as a function of the spectral properties of fully 
inhomogeneous magnetic fields. 

The present paper is then organized as follows. In Section 2 the initial condtions of 
the Einstein- Boltzmann hierarchy will be discussed for typical wavelengths larger than the 
Hubble radius prior to equality. In Section 3 the tight-coupling equations will be derived 
and discussed analytically. The numerical treatment of the system is presented in Section 4. 
Section 5 contains our concluding remarks. 

2 Adiabatic and non-adiabatic modes 

2.1 Zeroth-order tight-coupling approximation 

By combining the Thompson cross section (defined, in the following, as (Tt) with the ioniza- 
tion fraction and with the electron density rie, the differential optical depth is given by 
k' = ricXeCTTa/ao. Since prior to decoupling the photon mean free path is much smaller than 
the Hubble radius at the corresponding epoch, the inverse of k' can be used as an expansion 
parameter. This is the basic idea of the tight coupling approximation. The lowest order, 
the tight-coupling approximation implies 1/k,' (or, equivalently, that oo). In this 

limit the divergences of the peculiar velocities of photons (i.e. 6^) and baryons (i.e. ^b) are 
driven to a common value i.e. ~ 6'b = O^-i^. In the following, the consequences of the 
zeroth-order in the tight coupling approximation will be investigated in the presence of fully 
inhomogeneous magnetic fields. Consider, therefore, the case where large-scale magnetic 
fields are fully inhomogeneous and characterized, in Fourier space,by the following two-point 
function 

{B,{k,r)B^{p,r)) ^ ^ Phk) PB{k) 5('\k + p) , (2.1) 

where 

P^ik)-{si-^, P^ik) ^ A^{^^J ■ (2.2) 

In Eq. (2.2) kp is the "pivot" scale at which the spectra are normalized (see also, below, Eq. 
(2.35)). Since the magnetic fields are fully inhomogeneous the line element of the background 
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geometry will not be anisotropic ^ (as it happens, for instance, in magnetized Bianchi-type 
backgrounds) and it can be parametrized as 



ds' ^ a\r)[dr^ - dx']. (2.3) 

Consequently, the evolution equations of the background geometry will be solely determined 
by the total pressure and energy density of the system as 

^' = ^aVt, n^-n'^A7rGa''(pt+pt), p[ + 3n(pt + Pt) ^ 0, (2.4) 

where Ti = a' /a and the prime denotes a derivation with respect to the conformal time 
coordinate r. Denoting the longitudinal fluctuations of the metric (2.3) as^ 

d^Qoo = 2a^0, d^Qij = 2a^ip5ij, (2.5) 

and selecting the conformally Newtonian gauge [53, 54], the evolution equations of photons 
and baryons can be obtained in terms of the density contrasts (i.e. 5-y and 5b) and in terms 
of the divergence of the peculiar velocities (i.e. 9^ and 6'b): 

5; = AiP' - ^e^, 6'^ = - + «:'(^b - ^7), (2-6) 

5; = 3^^' - ^b, e'^ = -m^ - vV + ^J^^ + \^-^n'{e, - e^). (2.7) 

a Pb 3 pb 

^ — * — * 

The term J x B parametrizes the Lorentz force contribution as it arises in a conducting 
plasma which is, however, globally neutral (see, for instance, [45, 46]). In this situation, 
typical of magnctohydrodynamics (MHD), the possible contribution of electric fields is sup- 
pressed by inverse powers of the conductivity as it will be discussed later. By subtracting 
the two equations for the velocities introduced in (2.6) and (2.7), the combination {9^ — 6'b) 
obeys 

— » — * — * 

{e, - e^y + + {e, - e^) = + ne^ + ZiR^. (2.8) 

V 3pb/ 4 aVb 

The term at the right hand side of Eq. (2.8) can be viewed as a source term for the evolution 
of the difference between the divergences of the photon and baryon peculiar velocities. In 
spite of the magnitude of the source term at the right hand side of Eq. (2.8), {6-y — 9^) will 
quickly be driven to zero in the limit a^- ^ oo. Consequently, even if pre-decoupling magnetic 
fields are present, the lowest order tight coupling expansion still implies that 9^ —>■ 6'b as 
expected. 



^If the expansion is anisotropic because of a uniform magnetic field oriented along a specific direction, as 
in the celebrated Zeldovich models [43, 44], the angular power spectrum does not depend solely upon i but 
also upon m. These situations can be successfully constrained through WMAP data (see [32, 33, 34]) and 
references therein. 

"^In Eq. (2.5) and in what follows, 6s{...) denotes the scalar fluctuation of the corresponding quantity. 
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In the tightly coupled limit, the evolution equation of = ^7 = 6'b can then be swiftly 
obtained by summing up the equations for the velocity fields in such a way that the scattering 
terms cancel. The result of this procedure is, in Fourier space. 



:R^r^°"4(i + i?b)^" ' ' 4(i + i?b) 



+ ^TT^7b = -J7T—^^^l + k <t> + 77^^(^B - 4aB), (2.9) 



where i?b = (3/4)pb/P7- For numerical purposes (see Section 4) the late-time cosmological 
parameters will be fixed, for a spatially flat Universe, as ^ 
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2.47 X 10-^ a;b = 0.023, uj^ = 0.111, a;^ = a;b + uj^, (2.10) 



where cux — h'^^x and Q,/^ — 1 — fl^n', the present value of the Hubble parameter Hq will be 
fixed, for numerical estimates, to 73 in units of km/(secMpc). In terms of the fiducial set 
of parameters of Eq. (2.10) the baryon to photon ratio appearing in Eq. (2.9) becomes 

698 \ f i-^h \ f <^7 



^•^-17^^10:023^12.47x10-5; ' ^2-^^) 

where z + 1 — ao/ais the redshift (for numerical purposes the scale factor will be normalized 
in such a way that qq — 1; see below, Eq. (2.25)). In Eq. (2.9) we traded the MHD current J 
(appearing in the second equation of (2.7)) for a combination of the magnetic energy density 
and of the anisotropic stress. More specifically, it can be easily shown that 

3 V-[/xE] 1 6 PS . 

4 a% = V <7B - -V Ob, = — , ^Pb-^,- (2-12) 

Equation (2.12) holds since, in MHD, the total Ohmic current J is solenoidal (i.e. V - J = 0) 
and given by An J = V x [45, 46] (see also the appendix of [28]). In addition, the function 
(Tb is nothing but a parametrization of the anisotropic stress that appears in the spatial part 
of magnetic energy-momentum tensor, i.e. 

= -SpB + Ui, dimi = ip^ + p^) VVb, (2.13) 

where 

'""^'-f' "' = 4i?(^'^'-T*0- f^-"' 

In various context the force- free approximation is often employed [49. 50] (sec also [28] and 
references therein). Such an approximation amounts, in practice, to set {V x B) x B ^ 0. 

"^Thc values of the cosmological parameters introduced in Eq. (2.10) are compatible with the ones 
estimated from WMAP-3 [2, 3] in combination with the "Gold" sample of SNIa [47] consisting of 157 
supernovae (the furthest being at redshift z = 1.75). We are aware of the fact that WMAP-3 data alone 
seem to favour a slightly smaller value of Wm (i.e. 0.126). Moreover, WMAP-3 data may also have slightly 
different implications if combined with supernovae of the SNLS project [48]. The values given in Eq. (2.10) 
will just be used for a realistic numerical illustration of the methods developed in the present investigation. 
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In this limit, from Eq. (2.12) ctb In the present framework this assumption shall 

not be invoked. 

According to Eq. (2.12), the leading contribution of fully inhomogeneous magnetic fields 
to the scalar modes of the geometry is given by the magnetic energy density and by the 
divergence of the MHD Lorentz force. These two quantities, affect the perturbed Einstein 
equations and, in particular, the Hamiltonian constraint 

- m{H(p + = iT:Ga^[5pt + S,pb\, (2.15) 

and the total (scalar) anisotropic stress 

V\(l) -ij) = UnGa^iip, + p,)V^a, + {p^ + p^)VVb]. (2.16) 

Equations (2.15) and (2.16) follow, respectively, from the (00) and {i ^ j) components of the 
perturbed Einstein equations written in the conformally Newtonian coordinate system; the 
quantity SgPt is the (total) scalar fluctuation of the energy- density of the fluid sources. It is 
relevant to remark that, in Eq. (2.16), on top of the anisotropic stress of the magnetic field, 
the only source of the anisotropic stress of the fluid is provided by the neutrinos that arc 
collisionlcss after weak interactions fall out of thermal equilibrium for typical temperatures 
of the order of 0.8 MeV. In analogy with Eq. (2.13), the anisotropic stress of the neutrinos 
has been parametrized as 

= -SsPt + nj, did'Ui = (p, + p.)VV,. (2.17) 

The quantity 5sPt is the fluctuation of the total pressure. The last constraint on the dynamical 
evolution us derived from the (Oi) components of the Einstein equations and it is 



, ^ V • (^ X B) 



(2.18) 



The second term at the right hand side of Eq. (2.18) is the divergence of the Poynting vector. 
In MHD the Ohmic electric field is subleading and, in particular, from the MHD expression 
of the Ohm law 

- - (V X B)xB , , 

ExB^^ ^ . 2.19 

47r(7 

Since the Universe, prior to decoupling, is a very good conductor, the the ideal MHD limit 
[46] can be safely adopted in the first approximation (see also [51]); thus for cr — > oo the 
contribution of the Poynting vector vanishes. In any case, even if a would be finite but large, 
the second term is suppressed in comparison with the first term at the right hand side of 
Eq. (2.18); such a term is given by the sum of the divergence of the peculiar velocity of each 
fluid component weighted by the speciflc enthalphy (i.e. (pa + Pa) ioi a generic a- fluid): 

(pt + Pt)^t=E(Pa + Pa)^a, (2.20) 
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or, in more explicit terms, 

4 4 4 4 

{Pt + Pt)Ot = -PuOu + -p^O^ + PcOc + PhOh = -^pvOy + pc + -^7(1 + Rh)0^h- (2.21) 

The second equality in Eq. (2.21) follows from the zeroth order in the tight-coupling expan- 
sion and from the definition of the baryon to photon ratio (2.11). Every regular solution to 
lowest order in the tight coupling expansion has to satisfy the momentum constraint of Eq. 
(2.18). Such a requirement implies, in particular, a restriction on the form of the peculiar 
velocity fields for comoving wave-numbers larger than the Hubble radius (i.e. fcr -C 1). 

The evolution of the cold dark matter (CDM) component only feels the presence of the 
magnetic field through the modifications induced by Qb and u-q on of the two conformally 
Newtonian potentials, i.e. and ■0: 

5'^ = 3V^' - ^c, + = (2.22) 

A similar situation occurs for the evolution of the neutrino component where, however, 
the effect of the magnetic field enters the not only the Newtonian potentials but also the 
anisotropic stress that is constrained from Eq. (2.16). The relevant equations for the neu- 
trinos are then, in Fourier space, 

-\e,, el^^S,-k^a, + e^, al = ^e,-^kJ^,s. (2.23) 
o 4 io iU 

The quantity ^^,3 appearing in Eq. (2.23) is the octupole of the (perturbed) neutrino phase 
space distribution while a^, (i.e. the anisotropic stress of the neutrinos introduced in Eqs. 
(2.16) and (2.17)) is related to the quadrupole moment of J-"^ as ai, — ^1/2/2. The rationale 
for the inclusion of jFj^s in the fluid description is that neutrinos are collisionless for tem- 
perature smaller than 0.8 MeV when the Universe is still deep into the radiation-dominated 
regime. In this regime the neutrinos do not interact, and, consequently, they should be rather 
described by the collisionless Boltzmann hierarchy of the (perturbed) neutrino phase space 
distribution. To address this problem, the usual approach is to use an improved description 
where on top of the monopole and dipole of the (perturbed) neutrino phase space distribu- 
tion (connected, respectively, with and 9^), one also takes into account the quadrupole, 
i.e. and sometimes (as in the case of specific non-adiabatic modes) also the octupole. 

With all these ingredients in mind, the evolution of the longitudinal fluctuations of the ge- 
ometry is determined from the trace spatial components of the perturbed Einstein equations, 
and can be written as 

t/j" + n{(t>' + 2^^') + {7e + 2n')(t> + \v\(t> - V') = ^T^Ga\5,pt + 5pb), (2.24) 

where S^Pt is the fluctuation of the total pressure introduced in Eq. (2.17). In the present 
situation it will be mandatory to set initial conditions for the fluctuations of the geometry 
and for the various species of the plasma before radiation-matter equality but, also, after 
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neutrino decoupling. In the conformally flat case described by Eq. (2.3), the scale factor that 
interpolates between the radiation epoch and the matter-dominated epoch can be obtained 
by integrating Eqs. (2.4): 



a[T) = aeq 



where Qmo and Qro are evaluated at the present time and the scale factor is normalized in 
such a way that qq — 1. In Eq. (2.25) ti = (2/ifo)-\/aeq/^mo- In terms of Ti the equality 
time is 

Teq ^{V2- l)n = 119.07 f ^) Mpc, (2.26) 



0.134 

i.e. 2req ~ ti. In this framework the total optical depth from the present to the critical 
recombination epoch, i.e. 800 < z < 1200 can be approximated analytically, as discussed in 
[52]. By defining the redshift of decouphng as the one where the total optical depth is of 
order 1, i.e. K{zdec,0) — 1, we will have, approximately 



2;dec ^1139 



^ 0^0268 



0.04317 ' 0.6462 + 0.1125 In (nb/0.0431)' 



where consistently with Eq. (2.10), h = 0.73. From Eqs. (2.27) and (2.25) it follows that 
for 1100 < ^dec < 1139, 275 Mpc < Tdec < 285 Mpc. 



2.2 Magnetized adiabatic mode 

In the limit t/ti ^ 1, the solution of the whole system of governing equations corresponding 
to the adiabatic mode can be obtained as an expansion in powers of kr. When the given 
wavelength is larger than the Hubble radius before equality (i.e. fcr < 1 for r < Tcq )the 
solution of the coupled system introduced in the previous subsection can be written as^ 



5u — 5j — —20* 



Or — -^^(t>*i 



4 = (^b = 



2 



0.^ 



R^^B - R' 



2 

„ R'y 

+ 2 ^ 



0* — -Ry^B 



R 



ctb 



+ ^ - 2<7B 



( 2, \ ih 
l + -i?J + -^[4aB - Ru^b], 



R^ 



Iff 
6R. 



Rv 



■(V* 



(2.28) 
(2.29) 
(2.30) 
0*), (2.31) 



^In the following, all the perturbed quantities will be defined in Fourier space and will depend on the 
comoving wave-number k. 
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where the fractional contribution of photons to the radiation plasma, i.e. has been 
introduced and it is related to R„, i.e. the fractional contribution of massless neutrinos, as ^ 

R^^l-R^, R^ = ^—, r=-N^[ — ) =0.68l(^). (2.32) 



1 + r' 8 Vlly V 3 

In Eq. (2.31) (j)^{k) and ip^{k) denote the values of the Foiuicr modes of the longitudinal 
fluctuations of the geometry that are (approximately) constant in time. Equations (2.28), 
(2.29) and (2.31) deflne a solution that is adiabatic since the density contrasts in matter (i.e. 
5c and 5c) are exactly 3/4 of the density contrasts in radiation (i.e. 5^ and 5y). In the limit 
(vanishing magnetic energy density) and ctb — (vanishing magnetic anisotropic 
stress), Eqs. (2.28), (2.29) and (2.31) reproduce the standard adiabatic mode deep in the 
radiation dominated epoch [42, 57]. Prom the second equation reported in (2.31) it is clear 
that the neutrino anisotropic stress is partially counterbalanced by the magnetic anisotropic 
stress. This aspect is entirely due to the presence of the magnetic anisotropic stress and it 
has been verified in different situations involving also entropic fiuctuations [58]. 

In view of the forthcoming numerical applications, it is useful to characterize the adiabatic 
mode in fully gauge-invariant terms. In fact the quantity introduced so far are meaningful in 
the longitudinal gauge. The spectrum of curvature perturbations used as initial condition for 
the numerical evaluation of the brightness perturbations is often introduced [1, 2, 3] in terms 
of the curvature perturbations on comoving orthogonal hypersurfaces [55, 56], customarily 
denoted with TZ. In the longitudinal gauge we have that this quantity is simply given by 

In spite of the fact that TZ is here defined in terms of the variables of the longitudinal 
gauge, its value is invariant under infinitesimal coordinate transformations. In terms of TZ 
the longitudinal fluctuations of the geometry entering the solution of Eqs. (2.28), (2.29) and 
(2.31) become: 

107^, 2R^{AaB - Ru^b) 
15 + 4i?, (15 + 4i?,) ' 

^ _ 2(5 + 2i?.) _^?^ _ 

^ 15 + 4i?^ 15 + 4:RJ ^ ^ ' 

where 7?.* (A;) is the constant value of curvature perturbations prior to matter-radiation equal- 
ity. Within the notations followed in this paper, the spectrum of TZ^{k) is deflned as 

{n4k)nM)-j^'Pn{k)5^'Hk + p), Vn{k)^An(^^^ . (2.35) 



^In the present paper /i^O^ = 1.68 x 10-^ i.e., recalling Eq. (2.10), /i^ilro = h?^^ + /i^il^ = 4.15 x 10"^. 
According to Eq. (2.25) this impHes that 1 + Zoq — 3228. 9(wm/0. 134). Possible mass terms for the neutrinos 
in the mcV range will be neglected here but can be introduced with a modification of the neutrino Boltzmann 
hierarchy used to derive the relations reported in Eq. (2.23). 
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In Eq. (2.35) is the "pivot" scale (already introduced in Eq. (2.2)) at which the spectra 
are normalized. Typical choices of kp range from 0.05 Mpc"^ down to 0.002 Mpc"^ [1, 2, 3]. 

Equations (2.1) and (2.2) allow to determine the spectra of flB{k) and aB{k). In the 
present notations ^ 

(Qb(^)^1b(p1) = ^V^{k)S^'\k+p), {Mk)aB{p)) = ^Va{k)S^'\k + p), (2.36) 
where 

rUk) = V.{k) = S{e)nl,(^-^y\ (2.37) 

The power spectra of Eq. (2.37) follow from the definitions of Qb and ctb by using the 
two-point function of the magnetic field, i.e. Eq. (2.1). In fact, according to Eqs. (2.12) 
and (2.14), Qb and ctb are quadratic in the magnetic field intensity, and, therefore, lead to 
two mode-coupling integrals that have been estimated in in the nearly scale-invariant limit 
(i.e. £ < 1) where the functions J-{s) and Q{s) are determined to be 

^ 4(6-£)(27r)^- ^ 4(188 -4£^-66£)(27r)^- 

£(3-2£)r2(£/2)' ^ 3£(3-£)(2£+l)r2(£/2)- ^ ^ 

In the opposite limit (i.e. £ ^ 1) the mode-coupling integrals arising in the expressions of 
Qb and (7b are dominated by (small scale) diffusive effects ( see, for instance, [35, 38]) and 
the power spectra are sensitive to the specific scale of the diffusion (Alfven) damping. In 
Eq. (2.37) the quantities 



2 



^BL^^, PBL = |^, p,^a\r)p,{r) (2.39) 

have been introduced. Notice that Bl measures the magnetic field intensity smoothed over 
a comoving length-scale L. In fact, by using a Gaussian window function e~^^-^^/^ for each 
magnetic field intensity, the (real space) magnetic autocorrelation function at two coincident 
spatial points is 

Bl = (A(f,r)5^(f,r)) =^B(^)V(|)(27r)- (2.40) 

''In the present paper the calculations are consistently conducted within the parametrization introduced 

in Eqs. (2.1) and (2.2). Different authors use slightly different conventions which can be, however, easily 
translated in the present parametrization. For instance, the magnetic spectral index defined as n in [14, 
35, 38, 39] correspond, in the present notations, to e — 3. In [14, 35, 38, 39] the factor 2n'^/k^ appearing 
in Eq. (2.1) is omitted. These definitions are clearly conventional but we regard as preferable the notation 
adopted here since it agrees with the usual way of assigning the power spectrum of curvature perturbations 
in CMB physics. The other remark is that the magnetic fields considered here are not helical (see, for 
instance, [40] and references therein). This is not a limitation since it can be easily appreciated that the 
helical contribution does not affect the scalar modes which are the main subject of the present investigation. 
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where = 2ti/L. Equation (2.40) implies that the amphtude of the magnetic power 
spectrum appearing in Eqs. (2.1) and (2.2) can be directly characterized in terms of B\ as 

It should be mentioned, incidentally, that nearly-scale invariant magnetic energy spectra 
can be achieved both in some class of inflationary models [59] and in the case of pre-big 
bang models [60]. It is not our purpose here to endorse any specific model. The spirit of 
this investigation is more modest since we shall be content of understanding if large-scale 
magnetic fields can be consistently included in the picture of pre-decoupling physics and 
in the framework of the usual techniques commonly employed for the description of scalar 
CMB anisotropics. 

On top of the curvature fiuctuations on comoving orthogonal hypersurfaces, another class 
of gauge-invariant quantities is related to the density contrasts both total and partial (i.e. 
pertaining to each single species of the plasma). The total density contrast on uniform 
density hypersurfaces can be defined as a functional of the longitudinal degrees of freedom 

p't 

coincides, for wavelengths larger than the Hubble radius, with TZ. In fact, using the def- 
initions (2.33) and (2.42) into the Hamiltonian constraint of Eq. (2.15) it can be verified 
that 

Having introduced related to the total density contrast, it is also useful to introduce the 
following 4 variables 

C. = -V^ + |, C-y = -V^ + ^, ^^ = -^ + 1' Cb = -V^ + |, (2.44) 

which are interpreted as the density contrasts for each independent fiuid on uniform curvature 
hypersurfaces. 

The evolution for the density contrasts described by Eqs. (2.6), (2.7), (2.22) and (2.23) 
become, in terms of the variables introduced in (2.44), 

Equations (2.44) and (2.45) are a useful tool for the discussion of the non-adiabatic solutions 
of the system as it will now be explicitly shown. 
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2.3 Magnetized non-adiabatic modes 

The inclusion of fully inhomogeneous magnetic fields not only affects the adiabatic mode 
but also the other (non-adiabatic) modes whose physical features have been investigated 
in the conventional case (see, for instance, [58, 61, 62, 63, 64, 65] and references therein). 
The adiabatic solution obtained in Eqs. (2.28), (2.29) and (2.31) implies, in terms of the 
quantities defined in Eq. (2.44), that 

C = C7 = Cc = Cb = 7^* - ^^^B. (2.46) 

Equation (2.46) is not accidental. It implies, in fact, that the entropy fluctuations of the 
mixture of fluids are vanishing. The entropy fluctuations are indeed deflned as 

5ij = -3(Ci-Cj), (2.47) 

where the indices run over the four components of the fluid deflned previously. Equation 
(2.46) implies, therefore, that the the entropy fluctuations vanish for the adiabatic mode. 

The pre-decoupling system also admit solutions where entropy fluctuations do not vanish. 
These modes are called non-adiabatic. If there are entropy fluctuations, the perturbation of 
the total pressure, appearing in Eqs. (2.17) and (2.24), can be written as 

SsPt = cldpt + 5pnad (2.48) 

where 

Sp.^ = T^^Micl, - cy^ij, cl = ^. (2.49) 

^/i-Pt ij Pt 

In Eq. (2.49) Cgj and Cgj arc the sound speeds of each (generic) pair of fluids of the mixture. 
According to Eq. (2.48), the fluctuation in the total pressure may arise cither because 
of a inhomogeneity in the energy density or thanks to the chemical inhomogeneity of the 
plasma. Chemically inhomogeneous means, in this context, that the plasma is constituted by 
various fluids with equations of state that are different and that entail necessarily a spatial 
variation of the sound speed. As a consequence of Eqs. (2.48) and (2.49), the evolution of C 
can be directly obtained from the evolution of SsPt stemming from the flrst-order covariant 
conservation equation: 

— * — * 

S,p[ + m{5,p, + 5sPt) - 3V^'(a + Pi) + (Pt + Pi)Ot = (2.50) 

Neglecting the Ohmic electric held we do obtain the wanted evolution of C, i-e. 

'hi ^ Ti, /' 2 ^ \ j- 



C = ^5Pnad + —— - - UpB - ^. (2.51) 

Pt + Pt + Pt V 3/ 3 

By using Eqs. (2.51) and (2.44) together with the explicit form of the evolution equa- 
tions in the longitudinal gauge, the analytic form of the non-adiabatic contributions can be 
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obtained. In the following, as an example, the magnetized CDM-radiation mode and the 
magnetized neutrino entropy mode will be presented. For the case of the CDM-radiation 
mode the solution, in the limit r < ri and kr < 1 can be written as 







— 



5, = 4V^l(;^) -R 

3 

k^n , 



3 



r 

Tl 



3 
4 

k^T, 



R^^B, 



77b 



5b = 3V^i I 

(01 + + ^[Ru^B - 4(7b], 



kr f R^ 



4|^aB - 1^ 



B 



-H^aB + 



(V^i- 



Tl/ 



where 



15 + iR^ 
8(15 + 2i?,) 



5^ + — -R^r^B 

4 ' 



01 



15 - 4i?^ 
8(15 + 2R,) 



+ — i^^r^B 



(2.52) 



(2.53) 



In the case of this solution, the longitudinal fluctuations of the geometry are vanishing for 
r < Tl and the relevant entropy fluctuations, i.e. Scy and Scu are constant and proportional 
to S^:. In the limit Qb and ^ the solution (2.52) reduces to the conventional 
CDM-radiation mode. While in the standard case the octupole moment of the (perturbed) 
neutrino phase-space distribution vanishes, in the magnetized case it is proportional to kr 
and is negligible for length-scales larger than the Hubble radius prior to equality. In analogy 
with Eq. (2.35), the spectrum of entropy perturbations will be parametrized as 



2n 



{S,{k)S4p)) = —Vs{k)5^'\k + ^, Vs{k) = A, 



k ^ 

kr, 



(2.54) 



where kp is the same pivot scale introduced in Eq. (2.35) while Ug is the spectral index of 
the entropy fluctuation. In the case of the magnetized neutrino entropy mode, the solution 
can instead be presented as: 



= 00, 1p 



2 ' 



00 



SRyRi/ 



5^ — —200 + -S^Ri^ — QbRj, 



3{AR, + 15) 
5i, = -200 



Q 3 3 



—S^Rj 



QbRj, 



3 3 

5c = ~2^° ~ ^^6-^7, 



= ^00, 



^7b — 



k\ 



200 



4 ~ 

-S^Ri, + Ri,^Ib — 4(7b 



14 



200 - \s.{Ry + 1) + 4^aB - fiei^o 



00- 



(2.55) 



In the case of the neutrino-entropy mode ^ and it is proportional to a constant that 
has been denoted in Eq. (2.55) as <S*(A;); the spectrum of S*{k) can be defined in full analogy 
with Eq. (2.54). An interesting feature of the magnetized neutrino-entropy mode (common 
to the conventional case when magnetic fields are absent) is that the curvature perturbations 
on comoving orthogonal hypersurfaces (i.e. TZ defined in Eq. (2.33)) vanish for wavelengths 
larger than the Hubble radius before equality. This statement can be verified by using Eq. 
(2.55) into Eq. (2.33) with the result that 7^ = -^/; - 0/2 = for kr < 1. 



3 Higher-order tight couphng expansion 

After having taken the derivative with respect to the conformal time r of the first relation 
appearing in Eq. (2.6), Eq. (2.9) can be used to ehminate the baryon-photon velocity; the 
result will of this manipulation is: 



,2 



r// I '^^^b k' 



2 



■ HRh I, k. 



-[4aB-^^B], (3.1) 



3(i?b + 1) ' 

where, according to Eq. (2.12), the Lorentz force has been traded for a combination of 
the magnetic energy density and of the normalized anisotropic stress. Equation (3.1) stems 
directly from the zeroth order in the tight-coupling expansion. 

To discuss the polarization, we have to go (at least) to first-order in the tight coupling 
expansion. For this purpose, it is appropriate to introduce the evolution equations of the 
brightness perturbations of the /, Q and U Stokes parameters characterizing the radiation 
field. Since the Stokes parameters Q and U are not invariant under rotations about the axis 
of propagation the degree of polarization P — {Q'^ + U'^y^'^ is customarily introduced [6, 66]. 
The relevant brightness perturbations will then be denoted as Ai, Ap. This description, 
reproduces, to zeroth order in the tight coupling expansion, the fluid equations that have 
been presented in Section 2 to set initial conditions prior to equality. For instance, the photon 
density contrast and the divergence of the photon peculiar velocity are related, respectively, 
to the monopole and to the dipole of the brightness perturbation of the intensity field, i.e. 

= 4Aio and 6^ = 3k An. The evolution equations of the brightness perturbations can 
then be written, within the conventions set by Eqs. (2.3) and (2.5) as 



A^ + {ikii + k')Ai + ikjicf) ^ ijj' + k' 



Aio + iiV], - ^P2{iJ,)Sp 



(3.2) 



k' 



A'p + {ikii + k')Ap = ^[1 - P2{l^)]Sp, (3.3) 
v'^ + Hv^ + ik(f) + ^ [Qb - 4(7b] + ^ (^^b + 3i An) = 0. (3.4) 
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Equation (3.4) is nothing but the second relation obtained in Eq. (2.7) having introduced 
the quantity i/cfb = Oh- The source terms appearing in Eqs. (3.2) and (3.3) include a 
dependence on P2(a*) = (3/i^ — l)/2 ( Pe{n) denotes, in this framework, the i-th Legendre 
polynomial); ; n — k-n is simply the projection of the Fourier wave- number on the direction 
of the photon momentum. In Eqs. (3.2) and (3.3) the source term Sp is defined as 

Sp{k, r) = An{k, r) + Apo(A;, r) + Ap2(A;, r). (3.5) 

The evolution of the longitudinal fluctuations is dictated by Eq. (2.16) and by Eqs. (2.15)- 
(2.18). According to the definitions used to derive Eqs. (3.2), (3.3) and (3.4), the £-th. 
multipole of the brightness perturbations is defined as 

Pi{fi)Ai{k, fi, T)dfi = 2{-iYAu{k, r), Pe{fi)Ap{k, fi, T)dfi = 2(-^)^Ap^(fc, r). 

(3.6) 

Equations (3.2) and (3.3) con now be expanded in powers of Tc = |1/k'|. In particular we 
can write 

Au{k, t) = Au{k, t) + Tc5u{k, r), Ap^(A;, r) = Ap^(A;, r) + T^Sptik, r), 

= Vb(fc,T) +rAb(/c,r). (3.7) 

Equations (3.2) and (3.3) will be preliminarily phrased in terms of a Boltzmann hierarchy 
coupling together the different multipoles of the brightness perturbations. Then the obtained 
equations will be expanded according to (3.7). The result of the first step of this procedure 
is 

A;o + A;An = V', (3.8) 



A[^ + hzAn - ^Aio = ^</' - 



Aii + -^ 



(3.9) 



Ai2 + \kAi^ - \kAn = -^[9Ai2 - (Apo + Aps)], (3.10) 
iU 

Ai, + k'Au = 2^[^Ai(^-i) - + ^ > 2. (3.11) 

Equations (3.8), (3.9) and (3.10) are obtained, respectively, by multiplying each side of Eq. 

(3.2) by -Po(a*)) -Pi(a*) and P2{lj) and by integrating over jj, with the help of Eq. (3.6). The 
same calculation can be repeated in the case of Eq. (3.3); the result is: 

Ai>o + Api = - [Ap2 + Ai2 - Apo] , (3. 12) 
2 k 

Ai,i + -A;Ap2 - -Apo = -k'Api, (3.13) 

^P^ + ^A;Ap3 - \kAp^ = -^[9Ap2 - (Apo + A^)], (3.14) 
iO 

A'p, + n'Ap, = ^A_[£Ap(,_i) - (£ + l)Ap(,+i)], £ > 2. (3.15) 
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Having derived Eqs. (3.8)-(3.11) and Eqs. (3.12)-(3.15), the system can now be expanded 
to the wanted order in Tc. As anticipated, to zeroth order in the tight couphng expansion 
Eq. (3.1) can be reobtained. In fact, to lowest order in Tc, 

Mk, r) = -3iAn(A:, r), A;^, + kAn = V', (3-16) 

while, from Eqs. (3.12)-(3.15), 

Apo(A;,t) = 0, e>0, Ak(A;,t)=0, £ > 2. (3.17) 

Summing up Eq. (3.8) (multiplied by 3i) with Eq. (3.4) (multiplied by Rh) and eliminating, 
from the obtained relation. An in favour of Ajq (see Eq. (3.16)) the decoupled expression 
for the evolution of the zeroth-order monopole becomes, as expected, 

a;; + ^A, + , X. = (r + - + (3.18) 

'° i?b + l '° 3(i?b + l) l + i?b 3V 3(i?b + l) ^ ^ 

which is the same of Eq. (3.1) if wc recall that 4Aio = S^,. The presence of the magnetic 
field modifies the evolution of the zeroth order monopole and, consequently, according to Eq. 
(3.16), the zeroth order dipole. Such a modification also implies an effect on the polarization. 
In fact, to first-order in the tight coupling expansion, Apo 7^ and it is proportional to the 
zeroth order dipole. Recalling the notations of Eq. (3.7) and expanding Eqs. (3.8)-(3.11) 
and (3.12)-(3.15) to first order in Tc it can be verified that 

5 I 3 

Spo{k,r) = -Si2{k,T), Sp2{k,T) = -di2{k,T), Si2{k,T) = —kAii{k,T), (3.19) 
4 4 io 

which also implies 

Ai2(A;, r) ~ r^k, t) ^ ^ kT,An{k, r). (3.20) 

Since Tc grows very fast during decoupling, to have a better quantitative estimate of the 
effect the procedure outlined in [6] can be followed. By taking the first derivative of the 
definition of Sp and by using Eqs. (3.11) and (3.15) we can get, directly the following 
evolution equation for Sp 

•^p + ^^^'-^p^^^An. (3.21) 

Equation (3.21) can be used to estimate the source term in Eq. (3.3) and, therefore, the 
polarization can be determined as 

Ap(A;, Tdec) = -^(1 - A^')Aii(Tdec)A:(7dec:z:i, (3.22) 

where adcc — 70-80 Mpc is the width of the visibility function K'e~'^ which has a Gaussian 
form[67]. The quantity Xi appearing in Eq. (3.22) is 

X, = / e-^^dn / — e-TS«^ ~ 1.719. (3.23) 
Jo Jl X 
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Let us now perform an oversimplified estimate. Let us ignore neutrinos and magnetic 
fields and let us compute out of Eq. (3.2) the ordinary Sachs- Wolfe contribution assuming a 
constant adiabatic mode that propagates smoothly through the matter-radiation transition. 
The constant mode will be ~ —(3/5)7^* where 7?.* is the constant mode of curvature 
fluctuations. Then, from Eqs. (3.18) and (3.22)-(3.23) 

Al(A;, Tdec) ^ -Y COS (/cCsbTdec), Ap(A;, Tdec) = "0.17(1 - //^)((7decA;Csb)V'm Sin (/cCgbTdec)- 

(3.24) 

which implies that the cross-correlation between temperature and polarization oscillates as 
sin(2A;Csbrdec)- 

It is finally interesting to compute the expression of the characteristic damping scale of 
the fluctuations arising from the dispersion relations. From Eq. (3.4), by taking the Laplace 
transform, i.e. Vb(A;, r) — ^ v\j{k,uj) (and similarly for the brightness perturbations), 
we will have that the following relation holds to second order in the tight-coupling expansion 



VY,{k,uj) — — 3i Aii(A;,a;) 



l-i 



-1 12 



+ e»(0)+C»(QB,(7B). 



(3.25) 



In the standard derivation of the dispersion relations the metric fluctuations are neglected. 
This procedure is also justifled in the present case provided the spectrum of the magnetic 
energy density is not too steep in frequency. In practice this is the case for e sufficiently 
smaller than 1. In this case also the magnetic contribution can be neglected. In this situa- 
tion, we can Laplace transform also the other equations (i.e. for the monopole, dipole and 
quadrupole). Then the first-order tight coupling expressions can be used into Eqs. (3.9) and 
(3.25) to ehminate v^^. The resulting expression will then be 



-u;\Rk + l) + — + t— 
6 \k'\ 



16 
45 



= 0, 



(3.26) 



which implies that u) = uji + iu!2 where 



k 



3(l + i?b)' 



0U2 = 



kd 



kri 



dr 



2K'{f) 



16 



+ 



45(i?b(f) + l 3(l + i?b(r))2 



(3.27) 



In practice, for the typical set of parameters introduced around Eq. (2.10) A;^^ ~ 3 Mpc, 
which sets, grossly speaking, the typical scale of diffusion damping. 



4 Numerical analysis 

The evolution equations in the tight coupling approximation will now be integrated numeri- 
cally. The normalization of the numerical calculation is enforced by evaluating, analytically, 
the Sachs- Wolfe plateau and by deducing, for a given set of spectral indices of curvature and 
entropy perturbations, the amplitude of the power spectra at the pivot scale. Here is an 
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example of this strategy. The Sachs- Wolfe (SW) plateau can be estimated analytically from 
the evolution equation of TZ (or Q by using the technique of the transfer matrix appropri- 
ately generalized to the case where, on top of the adiabatic and non-adiabatic contributions 
the magnetic fields are consistently taken into account. The result of this calculation can be 
expressed, for £ < 30, as [7] 



— Zi{nr,£) + — R^njiL22{e,i) - AnAs2i{nrs,i) cos7^ 



9 



4 
25 



4 3 / — 

zo zo 

6 / — 

+—yAsRj^BL23{ns,e,£) cos 



(4.1) 



where the functions Zi, Z2 and Z^ 



Z3(n,e,e) 



.Vfeoy-y r(3-n)r^f + ^ 
4 V/cp/ r2(^2- f^r(^£ + | - f 

^\2e^(,.{koV' r(2-25)r(£ + 5 



r2(|-s)r(£ + 2-£) 



£ n 

2 4 



(4.2) 
(4.3) 

.(4.4) 



We recall the meaning of the various spectral indices appearing in Eqs. (4.2), (4.3) and (4.4). 
The spectral index e has been introduced in Eq. (2.2) (see also Eq. (2.37) for the magnetic 

energy spectrum); kg = Tq^ where tq is the present observation time. In Eqs. (4.2), (4.3) 
and (4.4), n stands either for Uf (i.e. the spectral index of the adiabatic mode defined in Eq. 
(2.35)) or for (i.e. the spectral index of the entropic mode defined in Eq. (2.54)) or for 
firs (i.e. the spectral index of the cross-correlation). In principle one could also introduce a 
different spectral index for each cross-correlation. However, following, for instance [65] we 
will assume that the spectrum of the cross-correlation is solely determined in terms of the 
spectrum of the autocorrelations, i.e., for instance, n^s = (n^ +ns)/2. In Eq. (4.1) 7^.5, 7^^ 
and 7sb are the correlation angles. From the request that the SW plateau is dominated by 
the adiabatic mode it is possible to constrain Qbl and, consequently, (see, indeed, Eq. 
(2.39)). The nature of the bound, however, depends on the parameters that enter Eq. (4.1). 
In the absence of magnetic and non-adiabatic contributions and for Eqs. (4.1) and Eq. (4.2) 
imply that for rir = 1 (Harrison-Zeldovich spectrum) £{£ + l)C^/27r = ^7^/25 and WMAP 
data [12] would demand that Atz = 2.65 x 10~^. 

If the SW plateau is determined by an adiabatic component supplemented by a (sub- 
leading) non-adiabatic contribution both correlated with the magnetic field intensity the 
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Figure 1: Bounds on the protogalactic field intensity as a function of the magnetic spectral 
index e for different values of the parameters defining the adiabatic contribution to the SW 
plateau. 

obtainable bound may not be so constraining (even well above the nG range) due to the 
proliferation of parameters. A possible strategy is therefore to fix the parameters of the 
adiabatic mode to the values determined by WMAP-3 and then explore the effect of a mag- 
netized contribution which is not correlated with the adiabatic mode. This implies, in Eq. 
(4.1) that As = and 7br = 7r/2. Under this assumption, in Fig. 1 the bounds on 
are illustrated. The nature of the constraint depends, in this case, both on the amplitude 
of the protogalactic field (at the present epoch and smoothed over a typical comoving scale 
L = 2T[/kL) and upon its spectral slope, i.e. e. In the case e < 0.5 the magnetic energy 
spectrum is nearly scale- invariant. In this case, diffusivity effects are negligible (see, for 
instance, [14, 15]). As already discussed, if e ^ 1 the diffusivity effects (both thermal and 
magnetic) dominate the mode-coupling integral that lead to the magnetic energy spectrum 
[14, 15]. In Fig. 1 (plot at the left) the magnetic field intensity should be below the differ- 
ent curves if the adiabatic contribution dominates the SW plateau. Different choices of the 
pivot scale kp and of the smoothing scale k^, are also illustrated. In Fig. 1 (plot at the left) 
the scalar spectral index is fixed to = 0.951 [3]. In the plot at the right the two curves 
corresponding, respectively, to rir = 0.8 and rir = 1 are reported. If e < 0.2 the bounds are 
comparatively less restrictive than in the case e ~ 0.9. The cause of this occurrence is that 
we are here just looking at the largest wavelengths of the problem. As it will become clear in 
a moment, intermediate scales will be more sensitive to the presence of fully inhomogeneous 
magnetic fields. 

According to Fig. 1 for a given value of the magnetic spectral index and of the scalar 
spectral index the amplitude of the magnetic field has to be sufficiently small not to affect 
the dominant adiabatic nature of the SW plateau. Therefore Fig. 1 (as well as other similar 
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plots) can be used to normalize the numerical calculations for the power spectra of the 
brightness perturbations, i.e. 



Let us then assume, for consistency with the cases reported in Fig. 1, that we are deahng 
with the situation where the magnetic field is not correlated with the adiabatic mode. It is 
then possible to choose a definite value of the magnetic spectral index (for instance e = 0.1) 
and a definite value of the adiabatic spectral index, i.e. rir (for instance rir = 0.951, in 
agreement with [3]). By using the SW plateau the normalization can be chosen in such a 
way the the adiabatic mode dominates over the magnetic contribution. In the mentioned 
case. Fig. 1 implies Bl < 1.14 x 10~^ G for a pivot scale kp — 0.002 Mpc~^. Since the relative 
weight of the power spectra given in Eqs. (2.35) and (2.37) is fixed, it is now possible to set 
initial conditions for the adiabatic mode according to Eqs. (2.28)-(2.31) and (2.34) deep in 
the radiation-dominated phase. The initial time of integration will be chosen as Ti = 10~^ri 
in the notations discussed in Eq. (2.25). According to Eq. (2.26), this choice implies that 

The power spectra of the brightness perturbations, i.e. Eq. (4.5), can be then computed 
by numerical integration. Clearly the calculation will depend upon the values of ujm, '^b, 
(Vc and R,^. We will simply fix these parameters to their fiducial values reported in Eqs. 
(2.10) (see also (2.11)) and we will take A^,^ = 3 in Eq. (2.32) determining, in this way 
the fractional contribution of the neutrinos to the radiation plasma. The first interesting 
exercise, for the present purposes, is reported in Fig. 2 where, in the left column, the power 
spectra of the brightness perturbations are illustrated for a wave- number A; = 0.1 Mpc~^; in 
the right column the power spectra of the same quantities are illustrated for k — 0.2 Mpc~^. 
Concerning the results reported in Fig. 2 different comments are in order: 

• for £ = 0.1 and rir = 0.951, the SW plateau imposes Bl < 1.14 x 10"^ G; from Fig. 2 
it follows that a magnetic field of only 30 nG (i.e. marginally incompatible with the 
SW bound) has a large effect on the brightness perturbations as it can be argued by 
comparing, in Fig. 2, the dashed curves (corresponding to 30 nG ) to the full curves 
which illustrate the case of vanishing magnetic fiels; 

• the situation where Bl > nG cannot be simply summarized by saying that the am- 
plitudes of the power spectra get larger since there is a combined effect which both 
increases the amplitudes and shifts slightly the phases of the oscillations; 

• from the qualitative point of view, it is still true that the intensity oscillates as a 
cosinus, the polarization as a sinus; 

• the phases of the corss-correlations are, comparatively, the most affected by the pres- 
ence of the magnetic field. 





(4.5) 
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Figure 2: The power spectra of the brightness perturbations for two typical wave- numbers. 
The values of the parameters are specified in the legends. The pivot scale is kp = 0.002 Mpc~^ 
and the smoothing scale is = Mpc~^ (see Fig. 1). 
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Figure 3: A detail of the cross-correlation (left plot). The autocorrelation of the intensity at 
Tdec as a function of e, i.e. the magnetic spectral index (right plot). 

The features arising in Figs. 1 and 2 can be easily illustrated for other values of e and for 
different choices of the pivot or smoothing scales. The general lesson that can be drawn is 
that the constraint derived only by looking at the SW plateau are only a necessary condition 
on the strength of the magnetic field. They are, however, not sufficient to exclude observable 
effects at smaller scales. This aspect is illustrated in the plot at the left in Fig. 3 which 
captures a detail of the cross-correlation. The case when = can be still distinguished 
from the case = 0.5 nG. Therefore, recalling that for the same choice of parameters the 
SW plateau implied that Bl < 11.4 nG, it is apparent that the intermediate scales lead to 
more stringent conditions even for nearly scale-invariant spectra of magnetic energy density. 
For the range of parameters of Fig. 2 we will have that Bl < 0.5 nG which is more stringent 
than the condition deduced from the SW plateau by, roughly, one order of magnitude. 

If e increases to higher values (but always with e < 0.5) by keeping fixed Bl (i.e. the 
strength of the magnetic field smoothed over a typical length scale L = 2TT/kL) the amplitude 
of the brightness perturbations gets larger in comparison with the case when the magnetic 
field is absent. This aspect is illustrated in the right plot of Fig. 2 where the logarithm (to 
base 10) of the intensity autocorrelation is evaluated at a fixed wave-number (and at t^^c) 
as a function of e. The full line (corresponding to a Bl = 10 nG) is progressively divergent 
from the dashed line (corresponding to Bl = 0) as e increases. 

In Fig. 4 the power spectra of the brightness perturbations are reported at Tdec and as 
a function of k. In the two plots at the top the autocorrelation of the intensity is reported 
for different values of Bl (left plot) and for different values of e at fixed Bl (right plot). In 
the two plots at the bottom the polarization power spectra are reported always at Tdec and 
for different values of Bl at fixed e. The position of the first peak of the autocorrelation 
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k = 0.002 MpG"',k|^ = 1 MpG-',n^ =0.951, e =0.1 



1^ = 30 nG 
., = 10 nG 




log(k Mpc) 



kp - 0.002 Mpc"' , k^- 1 Mpc"^ , = 0.951 , B^^ - 1 nG 




log(k Mpc) 



k =0.002 Mpc"', k^- 1 Mpc-',n^ = 0. 951,e-0.1 



kp -0.002 Mpc"', kj^- 1 Mpc"', = 0.951 , e - 0. 1 





log(k Mpc) 



log(k Mpc) 



Figure 4: The power spectra of the brightness perturbations at T^ec for the parameters 
reported in the legends. 
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of the intensity is, approximately, fej ~ 0.017 Mpc~^. The position of the first peak of the 
cross-correlation is, approximately, 3/4 of k^. From this consideration, again, we can obtain 
that Bl < 0.3 nG which is more constraining than the SW condition. 

Up to now the adiabatic mode has been considered in detail. We could easily add, 
however, non-adiabatic modes that are be partially correlated with the adiabatic mode. It 
is rather plausible, in this situation, that by adding new parameters, also the allowed value 
of the magnetic field may increase. Similar results can be achieved by deviating from the 
assumption that the magnetic field and the curvature perturbations are uncorrclatcd. This 
aspect can be understood already from the analytical form of the SW plateau (4.1). If there 
is no correlation between the magnetized contribution and the adiabatic contribution, i.e. 
Ibr = 7r/2, the SW plateau will be enhanced in comparison with the case when magnetic 
fields arc absent. The same situation arises when the two components are anti-corrclated 
(i.e. cos 75,. < 0). However, if the fluctuations are positively correlated (i.e. cos 75,. > 0) the 
cross-correlation adds negatively to the sum of the two autocorrelations of TZ and JIb so that 
the total result may be an overall reduction of the power with respect to the case 75^ = 77/2. 

5 Concluding remarks 

In the recent past, various attempts have been made to constrain large-scale magnetic flelds 
from CMB physics. While interesting results have been obtained for the vector and tensor 
modes of the geometry, the scalar modes have been, comparatively, less studied. In this 
study a systematic strategy to tackle this problem has been proposed and developed. The 
basic inspiration of the present approach can be summarized by stating that the effects of 
fully inhomogeneous magnetic flelds on CMB anisotropies is not univocal. In the standard 
lore the only relevant parameters to deflne the effects of large-scale magnetic flelds on CMB 
anisotropies are the ones related to the magnetic held itself. For instance, if the magnetic held 
is uniform the only parameter will be its intensity; if the magnetic fleld is fully inhomogeneous 
the only two parameters will be the spectral amplitude and the spectral slope. The point 
of view of the present investigation is that this way of thinking is ambiguous and, to some 
extent, not so productive. In fact, the spectral parameters of the magnetic field are certainly 
a necessary ingredient of the analysis. They are, however, not sufficient to conduct the 
calculation of the brightness perturbations. The information on the spectral parameters 
must be complemented by the specific pre-decoupling initial conditions. According to the 
results of this investigation, a sound procedure to be followed can be summarized as follows: 

• include the magnetic fields in the system of scalar modes before matter radiation 
equality; 

• solve the system in the long wave-length approximation prior to decoupling; 
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• use the obtained solutions for the numerical evaluation of the brightness perturbations, 
for instance, in the tight coupling approximation. 

The ambiguity mentioned above arises during the first step of the procedure we just outlined. 
In fact, the inclusion of magnetic fields prior to equality does not lead to a single solution 
but to different solutions. We will then have quasi-adiabatic and quasi-entropic modes where 
a magnetized contribution can be accommodated. 

In some sense, the present proposal is more conservative than previous attempts: the 
signature of magnetic fields on CMB physics is not directly accessible but it is mediated 
by the (supposedly rich) parameter space of the pre-decoupling initial conditions of CMB 
anisotropies. This approach is indeed common when mixtures of adiabatic and non-adiabatic 
modes are studied in the absence of magnetic fields. In the latter case the philosophy is to 
specify what kind of non-adiabatic mode we want to study and constrain. Here the situation 
is similar and we do hope, that, in the future, the considerations reported here will serve as 
an inspiration for analyzing the effects of fully inhomogeneous magnetic fields on the scalar 
CMB anisotropies. 

This analysis reported in this paper allowed to extend the tight-coupling approximation 
to the case when large-scale magnetic fields participate to the dynamics of the plasma. The 
evolution equations, in the tight coupling approximation, have then been integrated numeri- 
cally using, as normalization, the analytical estimate of the Sachs- Wolfe plateau modified by 
the presence of large-scale magnetic fields. The reported results allow not only to set con- 
straints but also to include a magnetized contribution in the current strategies of parameter 
extraction. 

From the combined analysis of the Sachs- Wolfe plateau and of the phases of Sakharov 
oscillations at intermediate scales, a fully inhomogeneous magnetic field uncorrelated with 
the adiabatic mode can be constrained to be smaller than 0.5 nC for a shghtly blue spectral 
slope e — 0.1 and for typical comoving wave-number of 1 Mpc~^. The conditions implied 
by the SW plateau ssem to be necessary but not sufficient since magnetic fields compatible 
with the SW bound may alter Sakharov oscillations at smaller length-scales. In the most 
general situation, correlation angles may proliferate: there may be correlations between the 
magnetized contribution and the adiabatic or even non-adiabatic modes. In this case the 
bounds on large-scale magnetic fields become much weaker and sizable magnetic fields may 
be allowed. 

The present analysis also shows that magnetic fields can affect the process of formation 
of CMB polarization. For instance, the usual approach to the Faraday rotation of CMB 
polarization is to assume that the polarization itself is only due to the adiabatic mode while 
the rotation of the polarization plane is only due to the magnetic field. This picture has been 
complemented by the observation that the magnetic fields affect the pre-decoupling initial 
conditions and, consequently, modify the evolution of the degree of polarization. 

The tight coupling expansion is rather useful to discuss the integration of brightness 
perturbations in some simple cases. Moreover, it is an essential ingredient of Boltzmann 
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codes since it is used, at early times, to avoid the integration of the (stiff) Euler equations. 
It will be interesting, for future studies, to analyze the space of magnetized initial conditions 
by solving directly the whole Boltzmann hierarchy. 
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